Superexchange-Driven Magnetoelectricity in Magnetic Vortices 
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We demonstrate that magnetic vortices in which spins are coupled to polar lattice distortions 
via superexchange exhibit an unusually large linear magnetoelectric response. We show that the 
periodic arrays of vortices formed by frustrated spins on Kagome lattices provide a realization of 
this concept; our ab initio calculations for such a model structure yield a magnetoelectric coefficient 
that is 30 times larger than that of prototypical single phase magnetoelectrics. Finally, we identify 
the design rules required to obtain such a response in a practical material. 
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The ability to control magnetism with electric fields, 
which can be realized through the interplay between spins 
and charges in solids, has an obvious technological ap- 
peal. The simplest form of such a control is the linear 
magnetoelectric effect, when an antiferromagnet placed 
in an electric field E becomes a magnetized with mag- 
netization M, while an applied magnetic field H induces 
an electric polarization P, proportional to the field: 

{Pi = a ijHj ri\ 
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Here onj is the magnetoelectric tensor and summation 
over repeated indices is implied. This magnetoelec- 
tric response requires simultaneous breaking of inver- 
sion and time-reversal symmetries, which defines the al- 
lowed magnetic symmetry classes and the non-zero com- 
ponents of the magnetoelectric tensor. While the phe- 
nomenology of the linear magnetoelectric effect is now 
well understood [H,0], the use of magnetoelectrics is ham- 
pered by rather low values of their magnetoelectric con- 
stants: for example in the prototypical magnetoelectric 
Cr2C>3, a (large) electric field of 1 x 10 6 V/cm induces a 
(tiny) magnetization of ~ 9 x 10 _5 /is per Cr ion. The 
search for materials with a much stronger response re- 
quires a deeper understanding of the microscopic mech- 
anisms of magnetoelectric coupling and new ideas about 
spin orders and crystal lattices that can conspire to pro- 
duce a large magnetoelectric effect. 

The aforementioned technological driver is likely also 
responsible for the renewal of interest in the related class 
of multiferroic materials which have simultaneous ferro- 
electric and (ferro) magnetic orders. Two recent develop- 
ments in this field are particularly relevant for the work 
we will present in this Letter. First, an early observa- 
tion that spiral magnetic order can lead to an electri- 
cal polarization^, has been confirmed repeatedly 0, Q 
and the list of such materials has been considerably en- 
larged. While spectacular non-linear magnetoelectric ef- 
fects, such as reorientation of electric polarization with 
a magnetic field, have been observed, the polarizations 



in such spiral magnets are small because the spin-lattice 
interaction is the weak spin-orbit-driven Dzyaloshinskii- 
Moriya interaction [(| At the same time, a new class 
of multiferroics has been identified in which the mag- 
netic ordering couples to the lattice through mechanisms 
of non-relativistic origin, in particular exchange stric- 
tion arising from superexchange0,0|. The stronger spin- 
lattice coupling leads to correspondingly larger magnet- 
ically induced ferroelectric polarizations, with polariza- 
tion values close to those of conventional ferroelectrics 
suggested. 

In this Letter, we show that these two concepts from 
the field of multiferroics - symmetry breaking in spiral 
magnets, and super-exchange mediated spin-lattice cou- 
pling - can be combined to yield materials with strong 
linear magnetoelectric response. 

We begin by considering the magnetoelectric response 
of a single spin vortex (Fig.QJ,). This can be viewed as a 
magnetic spiral rolled into a circle, and so we can use the 
results from spiral multiferroics to analyze its magneto- 
electric response. The magnetically induced ferroelectric 
polarization in spiral multiferroics is described by 

PocexQ , (2) 

where e is the axis around which the spins rotate and 
Q is the spiral wavevector[10|. In our context, this cou- 
pling induces an inhomogeneous electric polarization lo- 
cally oriented along the radial direction, so that the net 
polarization of the vortex is zero. A magnetic field ap- 
plied in the xy plane leads to a non-uniform rotation of 
spins in the vortex, which results in a nonzero net elec- 
tric polarization proportional to the magnetic field (see 
Fig. [TJd) . The spin vortex shown in Fig. [T^ has a diago- 
nal magnetoelectric tensor, with magnetization induced 
parallel to the applied electric field, while for the vortex 
shown in Fig. [TJ; an applied magnetic field induces a per- 
pendicular electric polarization and the magnetoelectric 
tensor is antisymmetric (see Fig. QJi). These conclusions 
are actually independent of the mechanism of magneto- 
electric coupling and generally follow from the fact that 
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the vortices shown in Figs.QJi and[T}; have respectively a 
monopole moment A oc ^ Q r a -S a and a toroidal moment 
Toc£ Q r Q xS a O. 
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FIG. 1: (a) A magnetic vortex carrying a pseudoscalar mo- 
ment. The thin solid arrows indicate the spin orientation, 
while the thick open arrows show the local polarization vec- 
tor; (b) A magnetic field applied to the vortex shown in (a) 
induces a net polarization along the field direction; (c) A mag- 
netic vortex carrying a toroidal moment; (d) A magnetic field 
applied to the vortex shown in (c) induces an electric polar- 
ization perpendicular to the field. 

Next we analyze the spin-lattice coupling resulting 
from the dependence of the Heisenberg superexchange in- 
teraction between spins on their relative positions. Con- 
sider a three-atom unit consisting of two magnetic tran- 
sition metal ions connected by a ligand such as oxygen 
that mediates superexchange (Fig. [2]). Due to charge dif- 
ferences, the cations and ligand shift in opposite direc- 
tions under application of an electric field. The exchange 
constant coupling the spins depends on the amplitude 
of the relative shifts through the changes in the metal- 
oxygen distance and the metal-oxygen-metal bond angle 
9. According to the Anderson-Kanamori-Goodenough 
rules[I3|, the exchange is antiferromagnetic (J > 0) for 
9 = 180° and ferromagnetic (J < 0) for 9 = 90°. Ex- 
periments varying A-site cation size in transition metal 
oxides have shown that the crossover from ferromagnetic 
to antiferromagnetic coupling is continuous [liij. There- 
fore, the total magnetization of the unit can be modified 
by applying an electric field; conversely changes in spin 
orientation will affect its electric dipole moment. 
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FIG. 2: (a) Two magnetic cations (solid circles) connected by 
a ligand (open circle); (b) Upon application of electric field 
the bond lengths and angle change resulting in a different 
relative alignment of spins Si and S2. 

Now we combine these two concepts to form a peri- 
odic array of magnetic vortices in which the magnetic 
moments are coupled through superexchange, and show 
that the combination leads to a large magnetoelectric re- 
sponse. The macroscopic magnetoelectric response of an 



array of magnetic vortices is proportional to the vortex 
density. Therefore we choose the smallest possible mag- 
netic vortex as our building block: a triangle of antifcrro- 
magnetically coupled spins, in which the angle between 
spins in the lowest-energy state is 120° (Fig. [3k.) . Using 
transition metal (TM) ions to provide the spins, and in- 
corporating oxygen ligands between them (Fig.[3h), leads 
to superexchange spin-spin interactions. Upon applica- 
tion of an electric held, the shifts of the oxygen anions 
relative to the positive TM ions induce changes in the 
Heisenberg exchange energy, changing the spin canting 
angles and resulting in a nonzero magnetization. The 
symmetry of the magnetoelectric response of the trian- 
gle is identical to that of the magnetic vortices of Fig. [TJ 
with the form of the in-plane magnetoelectric tensor con- 
strained by its C3V symmetry to 
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Here if is the angle between spins and the vectors directed 
from the triangle center to the corresponding vertex. In 
particular, <p = (Fig.[3h) leads to ev^ = — aoSij, so that 
for ao > the induced magnetization M is antiparallel 
to the electric held E; for ip = tr/2, M is perpendicular 
to E (Fig. Eh). 





FIG. 3: Magnetoelectric response of a single TM-O triangu- 
lar unit. For antiferromagnetic exchange coupling the angle 
between the spins Si, S2 and S3 (solid arrows) in the zero- 
field ground state is 120°, so that the net magnetization is 
zero. For <p — [panel (a)] all spins are oriented out horn 
the center of the triangle resulting in a nonzero pseudoscalar 
moment. Upon application of an electric field, E, the oxygen 
atoms (open circles) displace (open arrows) relative to man- 
ganese (solid circles) inducing a net magnetization through 
changes in the exchange coupling. The net magnetization M 
is then opposite in direction to E for if = 0, regardless of the 
orientation of E with respect to the spins. For f = \ [panel 
(b)] the spin triangle has a toroidal moment and the induced 
magnetization is perpendicular to the electric field. 

To transform the concepts outlined above into a model 
material with a three-dimensional periodic structure, we 
begin with planes of Mn atoms situated on the vertices 
of a Kagome lattice and assume that their spins form 
the 120° structure with zero wave vector (see Fig. 2J, as 
observed e.g. in iron jarosite KFe^(OH)fi(SO/i )o[l3|. At 
first glance, such a spin lattice would yield no magne- 
toelectric response because spins in the vortices formed 



at "up" and "down" triangles are oriented in opposite 
senses [if — and ip — n). However, when oxygen ions 
are positioned outside the "up" triangles and inside the 
"down" triangles, the sign of magnetoelectric coupling 
[otQ in Eq. [3] also alternates and the contributions of all 
triangles to have the same sign. This can be under- 
stood by comparing the magnetoelectric response of the 
S2 — Si and Si — S' 2 spin pairs and noting that, for fixed 
bond lengths and angle, only the scalar product of spins 
is important in Hciscnbcrg exchange. 

The two-dimensional plane shown in Fig. [4] has a sim- 
ilar structure to the MnO lay ers of the experimentally 
realized YMnC>3 structure 15| , which consists of a con- 
nected mesh of oxygen trigonal bipyramids with Mn 
atoms at their centers. Using this structure as mo- 
tivation, we extend our two-dimensional MnO planes 
to a three-dimensional periodic structure and introduce 
counter-ions (Ca and Al) in the voids of the lattice so that 
the correct charge balance is attained. To ensure that 
the sign of magnetoelectric response is the same for all 
layers, the neighboring MnO planes are rotated by 180° 
with respect to each other (Fig. [5|). This reverses the 
positioning of oxygen ions with respect to the "up" and 
"down" spin triangles in the next layer, which compen- 
sates the reversal of the spin direction that must result 
from the antiferromagnetic interlayer coupling provided 
by the 180° connections through the apical oxygen atoms. 
Our resulting 'KITPite' structure [22], with chemical for- 
mula CaAlMii307, correctly breaks / and T symmetry; 
in addition the apical oxygen ions between the MnO lay- 
ers are centers of combined IT symmetry. 




FIG. 4: The structure of one MnO plane. The Mn atoms 
(solid circles) are arranged on a Kagome lattice (dashed lines) 
with oxygen atoms (open circles) mediating the binding and 
superexchange. 

Finally, to assess the strength of the magnetoelectric 
response of the model material introduced in the previous 
section, we turn to first principles calculations employ- 
ing plane- wave density functional theory (DFT) as imple- 
mented in VASPflij]. We use PAW potentials for core- 




FIG. 5: Two layers of KITPite in the relaxed structure (space 
group Pmma). The Mn ions are at the center of the oxygen 
trigonal bipyramids (purple polyhedra). The correct charge 
balance is obtained through Ca (blue) and Al (green) counter 
ions situated in voids of the MnO mesh. 



valence separation 17], and include non-collinear mag- 
netism for the valence electrons. We approximate the 
exchange-correlation part of the Kohn-Sham potential 
using the rotationally invariant form of the LSDA+J7 
in the fully localized limit [l8| , with the Hubbard U 
applied only to the Mn d electrons (U = 5.5 eV and 
J = 0.5eV[19|). It is important to note that we delib- 
erately do not include spin-orbit coupling in our calcula- 
tions so as to ensure that our calculated magnetoelectric 
response arises entirely from the superexchange coupling. 
As a result, the pseudoscalar and toroidal spin arrange- 
ments (Figs. |3] (a) and (b)) are degenerate, and we are 
unable to determine whether a is diagonal or antisym- 
metric. 

We first relax the structure in the absence of an elec- 
tric field by optimizing the ionic coordinates to find the 
lowest-energy state with the constraint that the trigonal 
bipyramids are prevented from tilting. This constraint 
preserves the Kagome structure of the MnO planes which 
is essential for demonstrating the superexchange magne- 
toelectricity of this model. The resulting structure is 
shown in Fig. [5l Our calculated valence electronic struc- 
ture is, as expected, similar to that of YMn03 15]: The 
formal Mn charge is 3+, with four majority d electrons 
per Mn (d xy , d x 2_ y 2, d yz and d xz ), and an unoccupied 
minority channel providing a local moment of 4 /ib /Mn. 

Subsequently, we apply an electric field and calculate 
the linear response of the ions, which is sufficient for com- 
puting the linear magnetoelectric coefficient. The force 
on an ion upon application of an external electric field 
is determined by the Born effective charge tensor, Z*, 
through Fpi — Z^Ej, where F is the force, E is the 
applied electric field, /1 is an index denoting the ion, and 
i,j are spatial directions. The summation convention for 
repeated indices is once again employed. All elements of 
the Z* tensor are computed through derivatives of the 
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bulk polarization Z*^ — j^ 2 -, where P is calculated us- 
ing the Berry-phase approach [2^] for a small displace- 
ment in all degrees of freedom individually. The method 
is close to that employed by Iniguez[2l"j]. 

In order to obtain the first-order ionic response to the 
field, we use the force-constant matrix C^vj — jjf 1 -- 
Then, to linear order, the ionic displacements for a given 
force are found by inverting the force-constant matrix 
through SR„j = C~j jSFpi, so that the ionic response to 
an applied electric field is 5R V j = C~- -Z^-Efc. Finally, 
the total magnetization is calculated as a function of £7. , 
yielding the linear magnetoelectric coupling constant. 

Fig. [Fj] shows the calculated magnitude of the induced 
magnetization as a function of applied electric field. With 
a field of 1 x 10 6 V/cm, the ionic response leads to an 
average displacement of Mn atoms of 0.007A in the di- 
rection of the field and of O atoms of 0.005A against the 
field. The Born effective charges, Z*, have an in-plane 
average magnitude of +3.30e~ for Mn and — 2.26e~ for 
O. Using an equilibrium spin arrangement as shown in 
Fig. |H the spins S 2 and S 3 rotate by -0.1° and 0.1° 
respectively, leading to a magnetoelectric coupling coef- 
ficient of a = 1.10 x 1(T 5 JT -1 V -1 m -2 . Transformation 
to regularized CGS units yields «cgs = 4.15 x 10~ 3 . 
For a benchmark, we compare to the magnetoelectric 
response of Cr^Oa computed also within density func- 
tional theory[21|, a(Cr 2 03) = 1.3 x 10~ 4 in Gaussian 
units (in good agreement with the experimental value). 
Hence, our model system has a magnetoelectric coupling 
around 30 times larger than that of Cr2C>3. Since the 
spin-orbit coupling was not considered in this work, mag- 
netic anisotropics that determine the angle tp were ne- 
glected. We therefore cannot predict whether KITPitc 
would carry a pseudoscalar or a toroidal moment in the 
ground state. However, the strength of the magnetoelec- 
tric coupling resulting from Heisenberg superexchange 
[ao in Eq.([3])] is insensitive to tp. 




FIG. 6: Calculated magnetoelectric response of the model 
system using density functional theory and a linear fit. 



In conclusion, we have combined the concepts of mag- 
netically induced polarization in magnetic vortices with 
lattice-mediated coupling through the superexchange 
mechanism to demonstrate strong magnetoelectric cou- 
pling in geometrically frustrated antiferromagnets. We 
showed that such a mechanism can be studied using 
modern density-functional theory approaches with non- 
collinear spin density functionals and augmented with 
linear-response methods, and we explicitly calculated the 
magnetoelectric coupling of a model transition metal ox- 
ide. While the linear magnetoelectric response of our 
model compound is larger than that of any known single 
phase material, we anticipate that many further improve- 
ments are possible: In particular materials with larger 
polarizability through increased Z*s or reduced rigidity 
would be promising. We hope that this study will stimu- 
late the search for additional novel strongly- coupled mag- 
netoelectric materials. 
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